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Abstract 

We study the exact statistical mechanics of Lame solitons using a transfer matrix method. This 
requires a knowledge of the first forbidden band of the corresponding Schrodinger equation with the 
periodic Lame potential. Since the latter is a quasi-exactly solvable system, an analytical evaluation 
of the partition function can be done only for a few temperatures. We also study approximately the 
finite temperature thermodynamics using the ideal kink gas phenomenology. The zero-temperature 
"thermodynamics" of the soliton lattice solutions is also addressed. Moreover, in appropriate limits 
our results reduce to that of the sine-Gordon problem. 
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1. Introduction 



Systems described by one-dimensional nonlinear scalar fields and governed by nonlinear 
energy functionals are ubiquitous in physics. Frequently the governing dynamical equations 
of these fields admit large-amplitude, localized, robust, particle-like solutions - the solitary 
waves Q| . In anharmonic models of various condensed-matter systems j^, [J the presence 
of these nonlinear excitations is thermally controlled. It is therefore important to proceed 
to a statistical mechanics investigation of the system, considering the different excitations 
that are potentially present, as well as their interactions, since they contribute nontrivially 
to the thermodynamical properties. 

Over the years, extensive research has been carried out seeking the soliton solutions of the 
integrable sine-Gordon (sG) equation [J and its (non-integrable) counterpart, the double 
sine-Gordon equation (DsG) j^y,^. The statistical and thermodynamical properties for 
both the sG and DsG potentials were studied j?i Q ■ 

There have also been advances in the study of the hyperbolic analogues of these prob- 
lems, the sine-hyperbolic Gordon (shG) equation [9j, respectively the double sine-hyperbolic 
Gordon (DshG) equation (which is non-integrable). The statistical mechanics of double 
sinh-Gordon solitons was also addressed in jicl ]. 

However, the study of the elliptic analogues of the sG and DsG problems is less advanced. 
The elliptic generalization of the sG is the Lame equation fill and its single soliton and 
soliton lattice solutions have been worked out only recently [12J. One would surmise that 
the elliptic generalization of the DsG equation is the Associated Lame equation 13, 14, 15^ . 
Nevertheless, we found |l2j that this is not the case, and this generalization is still an open 
problem. 

There are many physical contexts in which the Lame equation arises, such as, e.g., bond- 
order and charge density wave systems 16 1 , nonlinear elasticity ^ | , phase slips in supercon- 



ductors and vortex oscillations in a Josephson junction [l^J, magnetoelastic interaction on 

appears in systems with 

a single periodicity, for instance a one- dimensional elastic solid comprising identical atoms 



with a certain strength of the potential |20L or Hartree equation for spinless fermions, or 
Hartree-Fock equation for spin-1/2 fermions [2l|. It also arises in field theory contexts such 
as sphaleron solutions of the Abelian Higgs model [3] , in connection with elliptic Calogero- 
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Moser systems [23], and in group theoretical aspects of algebraic potentials [24]. 

Our goal here is to address the problem of the statistical mechanics of thermally con- 
trolled excitations in a nonlinear system with a Lame potential. The presented results find 
applications in the various classes of physical systems enumerated above. Note, however, 
that the Lame equation is a quasi- exactly solvable problem, which means that one cannot 
obtain analytically the complete band spectrum of the corresponding Schrodinger equation, 
but only a part of the spectrum, and that only for certain values of the parameters of the 
potential, as it will be discussed below. Therefore, it is not possible to obtain analyti- 
cally the corresponding partition function for arbitrary temperatures, but only for a set of 
temperatures. 

The plan of the paper is as follows. In Sec. 2 we describe the soliton-like solutions of 
the Lame potential. For the sake of completeness, we review some previously obtained 



results 



12j and then present a new family of kink lattice solutions. We also address the 



zero-temperature (T = Oi^) "thermodynamics" of the kink crystals and discuss their sta- 
bility. The single topological kink solution is also described, as well as the asymptotic 
interaction energy between kinks. The sine-Gordon results are recovered in an appropriate 
limit. Section 3 describes an approximate method of treating the low-temperature statistical 
mechanics of this kink-bearing system, namely the ideal kink gas phenomenology. In Sec. 4 
we present the transfer operator formalism for an exact calculation of the grand-canonical 
partition function of the system, with some computational details relegated to the Appendix. 
Finally, conclusions and perspectives are presented in Sec. 5. 



2. Soliton Solutions of the Lame Potential 



The Lame potential |l3j, ll4j, 3] with a dimensionless parameter v is defined as 

V L {4>) = z/(^ + lH 2 sn 2 (0, k), (1) 

where v{y + 1) > for the situations of physical relevance. Here sn(0, k) is the sine- 
amplitude Jacobi elliptic function of real modulus k (0 ^ k ^ 1) and period 4K(k); K(k) 
denotes the complete elliptic integral of the first kind, see Refs. HQ- 

Lame potential is 

periodic with a period 2K(k), 



V L (<j> + 2K(k)) = V L (<f>) 



(2) 



and has one minimum V m i n = at <fi = and one maximum V max = v{y + l)k 2 at = K(k), 
see Fig. [TJ 

Consider a classical dimensionless scalar field <f>(x, t) with Vl(4>) as the interaction energy 
density. Its dynamics is described by the Euler-Lagrange equation 

&± - - - ^ n\ 

at 2 dx 2 ~ u ° d<\> [ } 

corresponding to a Hamiltonian density 

'1 



H(x, t) = A 



(4) 



Here V(x,t) = d<p/dt = <p t is the conjugate field momentum density, and <p x = d<fi/dx. In 
view of the periodicity of the Lame potential, the field is defined modulo 2K(k). The veloc- 
ity c, the frequency c^O; and the constant A (that has the dimensions [energy] -[time] 2 /[length]) 
are the parameters of the system defining the characteristic time (1/uq), length (c/ujq), and 
energy (Acu ) scales. In the following we shall use dimensionless quantities based on this 
scaling. 

The total energy 

E A [4>, V] = I H(x,t)dx (5) 



J —A 

and the topological charge (or the winding number) of the system 

are constants of the motion for fixed boundary conditions (2A is the length of the system, 
—A ^ x ^ A). In the "thermodynamic limit" A — ► oo, as seen below, the topological charge 
is simply the difference between the numbers of kinks and antikinks in the system. 
In the stationary case the equation for the field (f>(x) reduces simply to 

= OVl m 
dx 2 d<p ' U 

that can be integrated by quadratures, see below. In view of the covariant form of Eq. Q, 
the time-dependent solutions <p(x, t) are immediately obtained from 4>(x) by a "Lorentz 
boosting" to velocity v, i.e., x — > (1 — v 2 )^ 1 ^ 2 (x — vt). 

There are two families of kink lattice solutions of (J7j). From the perspective of the 
dynamical systems theory, the first family of non-topological kinks corresponds to closed 
periodic orbits in the {(f) , <p x ) phase plane, while the second family of topological kinks 
corresponds to open periodic orbits in the phase plane. These different kink lattices reduce 
to the same single soliton solution in the appropriate limit. 



2.1. First family of kink lattice solutions 



A first type of stationary solutions is given by 

±V2(x - x )= / y 8 

U(x ) V^W, k) - u(u + l)k 2 a 2 

where xq and a 2 are the two integration constants. We fix xo through the convention 
0(0) = K(k). The control parameter a 2 , with ^ a 2 < 1, determines the appropriate 
integration domain in Eq. (J8|). the profile and the properties of the field (f>(x). 

When extended to the whole real axis A — > oo, Eq. (JSJ) leads to a soliton lattice solution: 

sn 2 (0 jfe) = (l-^ 2 )-(l-a 2 )sn 2 (y, t) 

W > ] (1 - k 2 a 2 ) - k 2 (l - a 2 ) sn 2 (y, t) ' 1 } 



where t = y/JT— a 2 ) /(l — k 2 a 2 ) and y = ^Jlv{y + 1)A; 2 (1 — k 2 a 2 ) x. This represents an 
array of nontopological solitons, i.e., a periodic sequence of "small" (nontopological) kink- 
antikink pairs. In Fig.|2]we show a typical profile of this type of soliton lattice solution (i.e., 
the field as a function of x) for a fixed value of a 2 . The spatial size of one kink/antikink 
(which represents half of the spatial period of the lattice) is a function of the control param- 
eter a 2 

2L(a 2 ) = ^= 2K(t) (10) 
y/2is(is + l)k 2 (l-k 2 a 2 ) 

and the topological charge of such a "small" kink/antikink is, respectively, 

±[1 - sn- 1 (a 1 k)/K(k)}. 

The energy per kink/antikink of the lattice, as a function of a 2 , can be computed as 



E K L(a 2 ) 



hv{v + \)k 2 







\k 2 


-i) 



k 



^U(k 2 t 2 , 



1 - k 2 a 2 

where U(k 2 t 2 , t) denotes the complete elliptic integral of the third kind j25l. 12 



(11) 



2.1.1. Kink lattice "thermodynamics " at T = OK 

One can immediately compute the thermodynamic characteristics of the ground state 
of this classical kink-antikink crystal. The ground-state energy per kink is expressed by 
Eq. while the specific volume per kink is given by Eq. (|10|). Correspondingly, the 



5 



thermodynamic pressure is obtained as 



/ dE KL \ dE KL /d(a 2 ) , 1U22 , . 

p - = - [l>m J M = -3C2L57553J = -" ( " + 1)k a ' (12) 



re last equality we have used the properties of the derivatives of the elliptic 
2f|). Thus, the internal energy of the crystal increases with increasing volume 



(in deriving t 
functions [2i 

under these isothermal conditions. The equation of state can be obtained formally by 
eliminating the parameter a 2 between Eqs. (jlOj) and (|T2*|) . 

The chemical potential at zero temperature \Ikl is just the enthalpy per kink, 



Hkl = E KL + 2LP KL = 2 



2v(v + l)k 2 



k 2 a 2 



k 2 a 2 ^,, k' 2 



K(t)-—Tl(k 2 t 2 ,t) 



k 2 k 2 



(13) 



where the complementary modulus k' 2 = 1 — k 2 . Finally, the isothermal compressibility at 
T = OK is 

1 d(2L) 1 d(2L)/d(a 2 ) K(t) - (1/a 2 ) E(t) 
XKL ~ ~2L dP KL ~ ~2L dP KL /d(a 2 ) ~ My + V)k 2 {\ - a 2 )K{t) ' ( ' 

are, respectively, the first- and second-kind complete elliptic integrals of 



261 ] . One notices that Xkl can change sign upon variation of the param- 



Here K(t) and Em 
modulus t, see (25, 

eters k 2 and a 2 , the states with xkl < being "mechanically" unstable. The corresponding 
stability diagram of this kink lattice solution at T = OK in the plane of the parameters 
(k 2 , a 2 ) is presented in Fig. El 

2.1.2. The sine- Gordon limit 



As discussed in [12], the Lame potential (JTJ) reduces to the sine-Gordon (sG) potential in 
the limit 

k^O, z/(z/ + l) oo, with v{y + \)k 2 -> Q = finite. (15) 
In this limit the above nontopological kink lattice solution reduces to <P s g{ x )'- 



sin 



2 ((j) sG ) = 1 - (1 - a 2 )sn 2 (^/2Qx, VI - a 2 ) = dn 2 {^/2Qx, VI - a 2 ) , (16) 



where dn(y / 2Qx, a/1 — a 2 ) is the delta-amplitude Jacobi elliptic function 



kink/antikink width is 2L s c(a 2 ) = 2K(y/l — a 2 )/\/2Q, and the energy per kink/antikink 
E sGKL = 2^2Q[E(VT^) - (a 2 /2)K(VT^)] . (17) 



The pressure of the sine-Gordon kink lattice at T = OK is P s gkl = —Qa 2 , and the 



chemical potential \i s gkl = 2\^^[E(y/T^a?) — (a 2 /2)K{yl — a 2 )]. The isothermal com- 
pressibility 

K(VT^) - l/a 2 E(VT^) 
XsGKL 2Q(l-a*)K(VT^) ' 1 J 

is always negative, which means that this class of sG kink lattice solution is unstable (see 

the k -> limit in Fig. El and H) 



2.2. Second family of kink lattice solutions 

The second type of stationary solutions is given by 

±V2(x-x )= I y 19 

Ufa) VVlW, k) + u(u + l)k^a 2 

where xq and a 2 ^ are the two integration constants. Here again we fix xq through the 
convention 0(0) = K(k). One notices that this family of solutions can be obtained formally 
from the previous one, Eq. Q, by replacing a 2 by —a 2 , and thus: 

VX n 1 - sn 2 (y, t) 

where £ = ^(1 + Pa 2 )/(1 + a 2 ) and y = v / 2u(u + l)k 2 (l + a 2 ) x with a 2 ^ 0. This is 
a topological kink lattice solution, see Fig. 01 i.e., a succession of topological kinks (or 
antikinks), of kink/antikink width (half lattice period) 

2L{a 2 ) = 1 ; . 21 

/ 2z/(z/ + l)P(l + a 2 ) 



(22) 



The energy per kink/antikink of the lattice is found to be: 



Ekl = 2\ 



'2v{v + l)k 2 
1 + a 2 



2.2.1. T = OK "thermodynamics" 

Proceeding as in the previous case of the first family of kink lattice solutions, one obtains 
for the pressure of the ground state of the second type lattice solutions 

P KL = u{y + l)k 2 a 2 , (23) 



and for the corresponding chemical potential 

1 



fJ>KL 



l 2v(is + l)k 2 



k 2 



+ a' 



1 + a 2 

Finally, the isothermal compressibility 

k 2 K(t) + (l/a 2 )E(t) 



K(t)-^U(k 2 J) 



Xkl 



(24) 



(25) 



2v{v + l)k 2 (l + k 2 a 2 )K(t) 
which is always positive, and thus the kink lattice solution of the second type is stable at 
T = OK. 



2.2.2. The sine- Gordon limit 



In the sine-Gordon limit (|15j) . the above kink lattice solution reduces to 
sin 2 (0 sG ) = 1 - sn 2 (v/2Q(l + a 2 )a;, l/Vl + a 2 ) = cn 2 ( y/2Q(l + a 2 )x, 1/Vl + a 2 ) , (26) 



of spatial period 2L sG {a 2 ) = 2K(\j\J\ + a?)/y/2Q(l + o?\ where 



cn(y / 2<5(l + a 2 )x, + a 2 ) is the cosine-amplitude Jacobi elliptic function |25|, |26j. The 

corresponding energy per kink is 



E s n k T, — 2 



sGKL — *\\ — 9 

v 1 + a A 



2Q 



1 + a 2 )E(l/VTT^) - (a 2 /2)K(l/VTT^)\ , (27) 

the chemical potential is [l s gkl = 2y / 2Q(l + a 2 )E(l/Vl + a 2 ), and the pressure 
is PsGkl = Qa 2 . The positiveness of the isothermal compressibility XsGkl = 
E(l/yl + a 2 )/[2Qa 2 i^(l/v / l + o?)] confirms the thermodynamic stability of the sG kink 
lattice of the second type. These results are in agreement with those found previously in [7( 
for the sG potential. 

2.3. The single soliton solution 

For the sake of completeness we consider the single soliton solution as well. The corre- 
sponding single soliton solution can be obtained from either of the above soliton lattices in 
the limit of "infinite kink dilution" , i.e., when the period of the lattice becomes infinite, 
L, L — > oo, which corresponds to the control parameter a 2 — > 0. It is a topological, "large" 
kink/antikink of charge 

.dx = ±l, (28) 



2K(k) 



oo 

4>x< 
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given by 



sn(0, k) 



(29) 



'l + k' 2 sinh 2 (?/*) 

where y* = \j2v{y + 1) k 2 x and k' 2 = 1 — A; 2 . The characteristic "spatial extent" of such a 
soliton is ^ink = 2/ \/2v{y + l)fc 2 . In Fig. |S]we present the typical profile of a single kink. 



The energy (rest mass) of the static single kink/antikink is shown 

E K = ^2u(v + l) In ■ 



121 to be 



(30) 



In view of the covariant form of equation (J21), a kink moving with a velocity v has an energy 



E K (v) 



,2\- 



l ' 2 E K = {E 2 K + V 2 ) 1 



(31) 



where V = Ex v(l — v 2 )- 1 ' 2 is the "relativistic" momentum. 

One can also compute the energy of the interaction of one kink/antikink in the lattice 
with the rest of the lattice as U(a 2 ) = Ekl{(i 2 ) — Ek (for the first type kink lattice) or 
as If (a 2 ) = Ekl{& 2 ) — Ek (for the second type kink lattice). From these expressions one 
can deduce (through a Taylor expansion in the parameter a 2 ) the long-distance interaction 
energy between two kinks. 

Note that there is another, more general way to obtain the asymptotic interaction energy 



between two solitons - Manton's method [12 , 
kink (|2^j). that in our case reads (e.g., for x — > oo) 

2 



281 ] . Given the asymptotic shape of the single 



cxp 



y/2l/(v + l)k 2 X 



(32) 



one obtains for the asymptotic interaction energy between two solitons separated by a dis- 



tance D > l/y/2u(u + l)k 2 : 
U as (D) « 



8y/2v(v + l)k 2 

: 1-k 2 



exp (-y/2v(v + l)k 2 D) . 



(33) 



The positive sign corresponds to a repulsive interaction between solitons of the same topo- 
logical charge, while between a kink and an antikink the interaction is attractive. For the 
particular case of the first type of soliton lattice, in the limit of "high kink dilution" (a 2 1), 
one has 

16 



D = 2L(a 2 ) 



2z/(z/ + l)k 2 In 



a 2 (l - k 2 ) 



> l/^/2u(u+l)k 2 



(34) 
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and thus the asymptotic attractive interaction energy between a kink and an antikink in the 
lattice is 

U as (a 2 ) « -a 2 y/v(v + l)k*/2. (35) 

The same expression is obtained in the limit of "high kink dilution" for the second type of 
soliton lattice. Finally, there is no difficulty or surprise in obtaining the properties of the 
sG single soliton in the appropriate limit (|I5|L see also Q. 

Next, let us briefly consider the question of the behavior of the single kink solution of 
the Lame potential with respect to small perturbations of its shape. For small oscillations 
i/j(x) exp(iu;£) around the stationary kink solution (J2HJ) . the linearized dynamics equation 
reads 

d 2 



2v{y + \)k 2 [1 - 2(1 + A; 2 )sn 2 (0, k) + 3£; 2 sn 4 (0, k)] - cj 2 j ij){x) = . 



dx 2 

Consider the following change of variables: 

2 



(36) 



(1 - fc 2 )cosh(2y*) + (1 + k 2 ) ' 

and 



(37) 



ii=[(l-k 2 )z/2] s X (z), 8S Ji- — ig— , (38) 



4 8z/(z/ + l)P 



With these new variables, the stability equation (jSUjl becomes Heun's equation 29]: 
d 2 X (2s + 1 , 1/2 1/2 \ d X 



that is a su 
following 



dz 2 \ z ~^~z — l^~z — l/k 2 )dz^~ 

(g _ i/gXg + 3/2)z - (g + l)(2s - 1)(1 + k 2 )/{2k 2 ) 

+ ,(,-l)(,-l/P) (39) 

jject of current intensive research in mathematical physics j^. Further on, 
29j , one can make a new change of variables 

z = sn 2 (n, k) , f(u) = [sn(u, k)]^' 2 X (z(u)) , (40) 



and the Heun equation above reduces to the Schrodinger equation 

d 2 f(u) 

+ V(u)f(u) = hf(u) (41) 

with the "potential" V{u) = (15/4)k 2 sn 2 (u, A;) + (16s 2 -l)/[4 sn 2 (-u, k)] andafixed "energy" 
h = 9(1 + A; 2 )/4. 
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As explained in detail in, e.g., [3J], and also checked directly, these stability equations 
always admit the Goldstone uj = (s = 1/2) mode as a solution, thus restoring the transla- 
tional invariance of the original Hamiltonian (J3J) broken by the introduction of the kink into 
the system. In addition to the translation mode at uj = 0, although difficult to find exactly, 
there may exist additional bound states with nonzero frequency of the stability equations. 
These solutions represent internal oscillations of the kink (or shape modes), corresponding 
to harmonically varying kink shape around the kink center. They are strongly localized 
around the static kink. Besides these discrete bound state solutions, there exist continuum 
states (extended modes) corresponding to phonons (see below Sec. 3.2) that are scattered by 
the kink (asymptotically, i.e., far from kink's center, this scattering results simply in a phase 
shift as compared to the kink-free system). Note that the degrees of freedom of the bound 
states (internal modes of the kink) are obtained at the expense of the phonon field (whose 
local density is altered by the presence of the kink), i.e., the excitation of these internal 
modes can be seen as a 'capture' of phonons by the kink, see for a detailed discussion of 
this point. 

One could legitimately address the question of the linear stability in the case of the 
kink lattice solutions, too. However, the corresponding linearized evolution equation of the 
perturbation (although reminiscent of a Heun's equation) is more complicated and will not 
be discussed further. 



3. Finite-temperature thermodynamics: The ideal kink gas approximation 



3.1. The contribution of the kinks 



A first approach to the study of the finite-temperature T > K thermodynamics of 
the field J) can be discussed within the framework of the ideal kink gas approximation, 
see 0, 8,[n] for a more detailed description of the method. This approximation consists of 
assuming that, at sufficiently low temperatures, very few kink- like excitations are present, 
and thus the mean distance between them is very large as compared to their typical size ^kink- 
Therefore, these kinks and antikinks behave like independent, non-interacting particles, with 
an energy (when moving with velocity v) given by Eq. (jBH). Moreover, the interaction 
between individual kinks and the phonons (see below) is also neglected, i.e., kinks and 
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phonons contribute independently to the thermodynamic properties of the system. 

Let us concentrate first on the contribution of the kink gas. Suppose there are Nk 
kinks and antikinks in the system; then W\ = Nk — is the constant topological 
charge (winding number), Eq. (ffi|) . There is no restriction on the ordering of these kinks 
and antikinks along the system. Let N\ = Nk + = 2Nk — W\ be the total number of 
kink-like excitations in the system. Then the grand- canonical partition function is simply: 

oo 

S A (/3, //) = J2 eM^N A )Z A ((3) , (42) 

N A =0 

where fi is the chemical potential associated with the creation of a kink-like excitation in the 
system, T is the temperature, (3 = 1/(/cbT), and Z\(j3) is the canonical partition function, 

= w { /r dv exp + p2)i/2 ] y (43) 

(with h Planck's constant). Thus, finally 

S A (/3, fi) = exp [(AEkA/^K^Ek) exp(j3fi)] 

« exp {(2A/h)(27iE K /(3)^ 2 exp[(3(fi - E K )}} (44) 

with the app roximation holding for low temperatures (3 ^> 1 (Ki is the modified Bessel 
function |25j). This leads to a pressure 

P K = (2nE K /p 3 h 2 ) 1/2 ex P [P(fx - E K )\, (45) 

and to an average kink-density (i.e., number of kinks per unit length) 

n K = (OPk/O^t = {2ti E K / (3k 2 ) 1 / 2 exp[/% - E K )\. (46) 

The situation of physical relevance is that when there is no external constraint on the 
number of kink-like excitations in the system, i.e., the kink density is solely determined by 
the temperature; this amounts to setting \i = in the above equation, thus obtaining 

n K = (27rE K /(3h 2 ) 1 / 2 exp(-/3E K ) , (47) 

which shows that the kink density is exponentially small at low temperatures /3 ^> 1. 

The free energy density (i.e., the free energy per unit length of the system) of this ideal 
kink gas is then simply 

F K = -(l/P)n K , (48) 
12 



and the internal energy density 

U K = d{f3F K )/d(3 = [E K - 1/(2(3)} n K . (49) 
Finally, the specific heat per unit length is 

C K = k B [{(3E K - 1/2) 2 - 1/2] n K . (50) 



In analogy with the correlation functions for periodic potentials [a, |32j, e.g. 
(sin 0(0), sin <j>{x)) for the sine-Gordon case, an appropriate measure of the static field-field 
spatial correlation function is given by the following expression: 

C{x) = (sn(0(O), k) sn(0(x), k)) , (51) 

where (...) designates the usual thermal average (in agreement with the convention in Sec. II, 
4>(0) = K(k)). The choice of the function sn(0, k) for the study of the spatial correlations 
of the field is related to its "sensitivity" to the presence of a kink (antikink); indeed, if 
a kink (i.e., a "jump"' of 2K(k) in the field 0) is present on the segment of length x, this 
leads to a change in the sign of sn(0, k). The value of C(x) can be easily estimated using the 
following qualitative argument at low kink density the probability Prob(n, x) of finding 
n kinks in a segment of length x is simply a Poisson distribution of mean (n) = ukx, i.e., 

(flftx) 71 

Prob(n,x) = j — exp(-riKx) . (52) 

Correspondingly, sn(0(O), k) sn(<f)(x), k) = (— l) n sn(0(O), k) = (— l) n . Thus the correlation 
function is simply 

oo 

C(x) = J^(-l) n Prob(n, x) = exp{-2n K x) , (53) 

n=0 

which corresponds to a correlation length = 1/(2^) of the field <fi. 



3.2. Linear modes (phonons) 

Let us briefly turn to the description of the phonons, i.e., the small-amplitude harmonic 
vibrations in the Lame potential well. The equation of motion (JJJ) for the field <fi can be 
linearized around the minimum = of the potential Vl{4>), thus leading to the following 
dispersion relation for the phonons: 

u\ = q 2 + 2is(is + l)k 2 . (54) 
13 



All the frequencies are real for all q ^ 0. In the semiclassical approximation, the contribution 
of the phonons to the free energy per unit length of the system at a temperature T is thus 
given by 

^^"(isl^^^' (55) 

where Ax is the size of a spatial cell corresponding to the smallest-wavelength linear mode 
of the system, i.e., to the Debye cut-off in terms of frequency (usually, Ax corresponds 
to the unit-cell length of the corresponding discrete lattice of the system under study). 
Taking the continuum limit Ax — > leads to the well-known "ultraviolet catastrophe" . The 
contribution of the phonons to the internal energy density is then 

U ph = d{(3F ph )/d(3 = l/(/3Ax) , (56) 

while the specific heat per unit length is a constant, C p h = 1/Ax. No interaction of the 
phonons with the kinks is taken into account in the above thermodynamic expressions. 

The total pressure, free energy density, internal energy density, and specific heat per unit 
length of the classical field <fi is the sum of the independent contributions of the phonons and 
the ideal kink gas, respectively. In a more refined approach jjj , one can take into account the 
interaction between phonons and individual kinks as a 'trapping' of phonons by the kinks, 
i.e., the excitation of the internal bound modes of the kinks at the expense of the degrees 
of freedom of the phonon field. However, this results only in a transfer of energy between 
the linear modes (phonons) and nonlinear ones (kinks), i.e., it does not affect the global 
thermodynamic properties of the system inside the limits of the approximate ideal kink gas 
scheme. 

Of course, all the thermodynamic properties of the field in the low-temperature limit 
are dominated by the contribution of the phonons, since the kink density, Eq. (I4*T[ ), is 
exponentially small at low temperatures. However, as discussed, for example in [3J, the 
density of kinks may be important for some transport properties that are insensitive to 
phonons, e.g. dc conductivity in charge- density- wave systems and transport coefficients in 



general 



3311 . which are related to the correlation functions of the field. 
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4. Finite-temperature thermodynamics: The transfer operator formalism 



The above approximate theory becomes inappropriate at high temperatures /3 < 1, when 
one can no longer neglect the interaction between the various nonlinear and linear excitations 
present in the system. There exists an exact general formal method for computing the 
grand-canonical partition function, by mapping the problem onto the spectral problem of a 



31 



Schrodinger operator. Let us briefly apply this so-called transfer operator formalism 
I32I to the case of the Lame potential. 

The canonical partition function for the field 4>(x) can be written as a functional integral 
over the conjugated fields 0(x, t) and V(x, t) = dcp/dt: 

Z A ((3) = Jw jv4> exp{-/3£ A [0, V}} . (57) 

One can define the path integral by dividing the system into M cells of length As = 2A/M 
and replacing the continuous fields <p(x, t) and V(x, t) by two discrete sets of (M + 1) field 
variables = <p(x u t)} and {P^t) = V(x l , t)}, x { = -L + 2iL/M, i = 0, ...,M. At the 

end, one would consider the continuum limit M — > 00. We shall consider fixed boundary 
conditions 

Mt) = 0, Mt) = 2K(k)W A ; V Q (t) = V M (t) = , (58) 

that assure the constancy of both the total energy E\ and the topological charge W\. 
Moreover, we impose that W\ is an integer, i.e., in the "thermodynamic limit" A — > 00 
the system supports topological kinks/ ant ikinks, and W\ represents simply the difference 
between the numbers of kinks and antikinks. 

The partition function factorizes into a kinetic energy part Z^(j3) (determined by the 
momentum variables {Pi}): 

l£) ' (59) 

and a configurational (potential energy) part Z^(P). The evaluation of the latter can be 
carried out using the formalism of the transfer operator j?l 31, 0], that allows an exact 
mapping of this problem onto the problem of finding the eigenvalues of an integral operator. 
In the continuum limit M 3> 1 (Ax A), this problem can be further simplified to finding 
the energy eigenvalues of the following Schrodinger equation for the Lame potential, with a 
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temperature-dependent "mass" : 

1 d 



2 ^ 2 |^ + W)-^]$n(0) = O, (60) 

where $ n (0) are the Lame functions jlll. Ilfil Is^. More precisely, 

/? A \ M//2 

Z*M = f-^J Eexp(-2A^ n )$ n (0)$;(2ir(A;)^ A ), (61) 

where (...)* denotes the complex conjugate, and n labels the "quantum states". 

In view of the periodicity of the Lame potential, Eq. (J2J), the eigenvalues of the Schrodinger 
equation ()6()j) lie in allowed bands, separated by forbidden bands. In each such allowed band 
the energy varies continuously with the wavenumber 

k = , where - 1/2 < q < 1/2 (62) 

K[k) 

(k is restricted to the first Brillouin zone). In the "thermodynamic limit" A — > oo, and using 
Bloch's theorem for the eigenfunctions, one can simplify further the expression (|61|) to: 

/27rAx\ M/2 f 1/3 

-1/2 

where $«j(0) is the eigenfunction in Eq. (jHOJ) corresponding to the eigenvalue £i(q) that lies 
within the first allowed band. 

The total canonical partition function in the continuum limit M ^> 1 is thus given by: 

/2tt\ m r 1/2 

Z A ({3) = ^— J J ^ dqexp[-2A/3e 1 (q)-2niqW A ]\$ q (Q)\ 2 . (64) 

One cannot pick out directly the most significant contribution to the integral above as the 
term of largest magnitude, because the phase 2nqW A also plays a role. The way to avoid 
this problem is to go to the grand canonical ensemble, whose partition function is: 

oo 

H A (/3,//) = exp(2A/3P) = Yl Z A (P)exp(fi(3W A ) = 

Wa=—oo 

r>^\M -1/2 



r27rAx\ M/2 r 1/2 

Zt((3) = (^p-) J i 2 ^exp[-2A/3 £l (g) - 2mqW A }\%(0)\ 2 , (63) 



P 



rfgexp[-2A/5e 1 (g)]|$ g (0)| 2 £ exp[-(27rig - nP)W A ] ■ 

1 / 2 W A =-oo 

(65) 



Here /i is the chemical potential associated with the creation of one topological charge in 
the system, and P is the thermodynamic pressure. Recall that W A represents the difference 
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between the number of topological kinks and antikinks in the system, and that is why it 
runs between — oo and +00. Allowing the chemical potential fi to take imaginary values 



H/3 = 2m\ , (66) 

one can immediately perform the summation over Wa, given that 

00 

exp[-2m(q - X)W A ] = 5(q - A) . (67) 

Wa=— 00 

Thus 

/2-7T \ M 

SaOM) = f j J exp[-2A/3e 1 (A)]|$ A (0)| 2 , (68) 
from which one obtains the pressure 

™-M^)-« m -M*w)-<-%)- (69) 

Note that the first term in the r.h.s., ^S^ln ( 27T ^ X j (corresponding to Vl = 0), is nothing 
else but the contribution of the classical phonon field, see also ^, yl|. 

We must now analytically continue the chemical potential back onto the real a xis, which is 
the situation of physical relevance. This problem has been analyzed in detail in [35] for gen- 
eral periodic symmetric potentials, and it was shown that under rather general circumstances 
the first allowed energy band of complex wavenumber maps onto the first forbidden energy 
band of real wavenumber. One is thus led finally to the important conclusion that in order to 
obtain the thermodynamics of the system, Eq. ([69)1 . one has to compute the lowest allowed 
and forbidden energy bands corresponding to the associated Schrodinger equation ()60|) with 
the Lame potential. 

This is a nontrivial spectral problem, and some details are presented in the Appendix. 
No general results are known for arbitrary values of (3 and v. However, for the case when 
2f3 2 v{v + 1) is of the form £(£ + 1), with £ an integer that is either strictly positive, or £ ^ —2, 
the energy bands are known under an implicit form that allows for combined analytical and 
numerical calculations. This means that, for a fixed value of u, the thermodynamic properties 
of the field can be computed for a discrete set of values of the temperature given by: 

"=VSS' ^^{0,-1}. (70) 



17 



As an example, we give below the analytical expression for the set of parameters with 
2f3 2 v{v + 1) = 2 (i.e., t = 1). In this case, the expression (p3§)) for the pressure of the system 
becomes: 

My/u(u + l) , ( v ! I ; 



2A 



In 



h 



1 



cn 2 (<5, k) 



(71) 



Here k' = y/l — k 2 is the complementary modulus of the elliptic functions, and the value of 
the parameter S follows from the implicit equation: 

5 K(k) ' 



7T 



F(7, k) - £7(7, k) + tan7 y 1 — Psin 7 



(72) 



2ti^/v(v + 1) 2K(k') 
where -£(7, k) and £7(7, k) are incomplete elliptic integrals of the first and second kind 
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261 ] . respectively, and 



sury = sn(<5, k) 



(73) 



From it one can compute, in principle, the thermodynamic properties of the system at this 
fixed temperature, e.g., the density of kinks, 

-1 



n K 



\dfi /T 



k' 2 ^v(v+l) sn(5, fc)dn(<&, k) 



1 



K(k) dn"(£, k) 
_2K(k') f ~!T cn 2 (5, k) 



(74) 



2tt C n 3 (5, k) 

One has now to consider the \x — > limit in the above expressions, that corresponds to the 
physically relevant situation of thermally activated excitations of the field, i.e., no external 
constraint on the topological charge of the system. 

To complete the discussion of statistical mechanics, we need to compute the field-field 
correlation function (|51|) whose Fourier transform is related to the static structure factor 
[3 Q, Q] • The corresponding correlation length is a measure of average distance between 
kink excitations. As discussed in 0, ^] for general periodic p_otentials, C(x) is obtained from 



the knowledge of the lowest Lame band wavefunctions [16l |24j , and in the asymptotic limit 
\x\ ^> 1 it is C(x) ~ exp(—2riK\x\), in agreement with the expression found in the low 
kink density limit, Sec. III. Its general expression, however, cannot be obtained in a closed 
analytical form. Similarly, the correlation function involving field fluctuations P, [3] for 
the Lame case remains an open question. Finally, as a general remark [3], the results of 
the transfer matrix formalism reduce in the low-temperature limit, through a WKB-type of 
approximation, to the ideal kink gas results in Sec. III. 
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5. Conclusions 



We have obtained two types of kink lattice solutions of the Lame equation and their 
common limit of the single-kink solution. We studied the T = OK "thermodynamics" of 
these kink crystals. Using the ideal kink gas formalism, we studied first the approximate 
low-temperature thermodynamics, which takes into account separately the contribution of 
independent kinks and phonons, while the interactions between them are not properly ac- 
counted for. Then we invoked the transfer integral approach to calculate exactly the par- 



tition function. This maps the statistical mechanics prob 



em onto the spectral problem of 



0. 



Unlike the sine-Gordon equa- 



the Schrodinger equation for the Lame potential, see also 
tion, Lame equation is not exactly solvable but quasi-exactly solvable: the band structure 
of the corresponding Schrodinger equation is not known analytically for all the values of 
the parameter I of t he p otential. The band structure of the Lame equation is known in 
explicit form for I — 1, and in implicit form for higher integer values of £, see 

the Appendix. Due to this limitation, we have been able to obtain closed form expressions 
for thermodynamic quantities only for a set of temperatures. Similar constraints hold in the 
case of another periodic system, namely the double sine-Gordon equation 0. Although we 
have not been able to present analytic results for all temperatures, our approach provides 
insight into this system and other related periodic quasi-exact solvable systems. 
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APPENDIX A: SPECTRAL PROPERTIES OF THE LAME POTENTIAL 

Let us consider the dimensionless Schrodinger equation 

d 2 



+ v L (<p) 



^(0) = • (Al) 
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The Lame potential with v = £ a strictly positive integer is well-known to be exactly solvable, 
see, e.g., jlll l34 j . The energy e and the wavenumber k are expressed parametrically as 



i 1 



n=l 



srrav 



E 

n=l 



cna n dna n 



sna 



n 



and 



n=l 



2K(k) 



respectively, where Z(u, k) is the Jacobi's zeta function (25j,|26j , 

Z(u, k) = E (sin-^sn^, k)), k) - (E(k)/K(k)) 



u 



(A2) 



(A3) 



(A4) 



The parameters ax, ■■■ ae are subject to (£ — 1) independent constraints 

D 

sna p cna p dna p + sna n cna n dna n 



E 

p=i 



sn z a„ — sn z a r 



0, p^n, n = !,...,(£-!). (A5) 



The allowed energy bands correspond to a real value of the wavenumber k, i.e., to the 
condition 



Re 



^Z(a m k) 



n=l 



o, 



(A6) 



and there are (£ + 1) allowed bands (£ finite bands followed by a continuum band), separated 
by £ forbidden bands. 

Closed analytical results can be obtained in the £ = 1 case. In particular, for the lowest 



allowed band one obtains the parametric equations 



24, 



34j for the energy 



e x = 1 - k' 2 cn 2 (r] + K(k'),k') , k'^Sx^l, ^ r] ^ K(k') 



(A7) 



and the corresponding wavenumber (for half of the first Brillouin zone, the other half being 
obtained by symmetry) 



k = Z(r), k') — 



7T 



(ex-k 2 )(l-ex) 



(A8) 



2K(k) \~ K(k') J V l + k 2 -ex 
Here k' = y/1 — k 2 is the complementary modulus of the elliptic functions. With a change 



of parametrization and using the properties of the elliptic functions [25 



26 



can be reduced to the following (more familiar to physicists) expressions 16l l2lj|: 



, these formulae 



El = l _ k' 2 cn 2 {p, k'), k 2 ^ ex ^ 1 , —K(k') ^p^ K(k') , 



(A9) 
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and for the wavenumber in the first Brillouin zone 



i=z M + mk*)- (A10) 



Note that (as demonstrated in a more general context in |35j]) the results for the first 
forbidden band can be obtained formally from those concerning the first allowed band simply 
by replacing p in (|A9J) and (|A10J) by ip. One obtains thus for the energy levels of the first 
forbidden band: 

k' 2 

1 cn 2 (p, k) ' ^ 

where the parameter p (-K(k') ^ p ^ K(k')) corresponds to a complex wavenumber given 

by: 

k = 



F(j, k) — -E(7, k) + tan7 y 1 — A; 2 sin 7 



+ t 2K(k)K(k>)> (A12) 



with 



sury = sn(p, fc) . (A13) 



One realizes, however, that with increasing i the complexity of the implicit equations 
(|A2j) . (|XH1),(|MJ), and (jMI) (that one has to solve in order to obtain the dispersion relations 
for the allowed bands) increases, up to becoming prohibitive even for present-day powerful 
computers. 



Fortunately, a recent paper by Maier 



3611 offers an alternative to the above- described 



Hermite-Halperin method of band computation. It is shown, using the so-called Hermite- 
Krishever Ansatz, that the dispersion relations for any integer I can be expressed in terms of 
the t = 1 relations. Moreover, the effective construction of these dispersion relations, whose 
technical details will not be given here, makes use of relatively simple polynomials and 
rational functions. This represents an enormous gain in the efficacy of symbolic calculation 
of the dispersion relations. 

To our knowledge, there is no information yet on the spectrum of the Lame potential 
when v is not a positive integer. In these cases one has to obtain the results strictly from 
numerical evaluations. 
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FIGURE CAPTIONS 



FIGURE 1: Typical profile of the Lame potential Vl(4>) over a period [0, 2K(k)]. 

FIGURE 2: Typical profile of a kink lattice solution of the first type (an array of nontopo- 
logical kink-antikink pairs) over a spatial period 4L of the lattice (for a fixed value of the 
control parameter a 2 , < a 2 < 1.) 

FIGURE 3: The stability diagram of the kink lattice solution of the first type, at T = OK, 
in the plane of the parameters k 2 (of the Lame potential) and a 2 (the control parameter). 

FIGURE 4: Typical profile of a kink lattice solution of the second type (an array of 
topological kinks/antikinks) over a spatial period 4L of the lattice (for a fixed value of the 
control parameter a 2 > 0). Note that we have plotted \<ft(x)\ modulo 2K(k); there are no 
"cusps" in the field 4>(x). 

FIGURE 5: The profile of the topological single kink solution. 
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FIGURE 2 
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FIGURE 4 
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FIG. 5: 
FIGURE 5 
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